A giant overshoot effect in the Janssen granular column 
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We present new experimental results on the mechanical behavior of static granular assemblies 
confined in a vertical column. Our measurements confirm, for the first time, the universal Janssen's 
scaling for the stress saturation curve. We show consequently, in the context of isotropic elasticity, 

a relation between the Poisson ratio and granular packing fraction. Moreover, using a systematic 
study of the overshoot effect created by a top mass equal to the saturation mass, we show behaviors 
reproduced qualitatively by isotropic clastic materials but in the case of a granular assembly of 
a spectacular amplitude. These experimental results are strong tests for any theory of granular 
matter. 



PACS numbers: 45.05.+x,45.70.Cc,46.25.-y 

The mechanical status of granular matter is presently 
one of the most open and debated question [1]. This 
state of matter exhibits many unusual mechanical and 
rhcological properties such as stress induced organiza- 
tion at the microscopic [2] or at the mesoscopic [3] level. 
This issue sets fundamental questions relevant to the un- 
derstanding of many other systems exhibiting jamming 
such as dense colloids or more generally soft glassy ma- 
terials [4, 5]. For practical applications, the quasistatic 
rheology of granular assemblies is described using a phe- 
nomenological approach, based on an elasto-plastic mod- 
elling of stress-strain relations [6]. So far, there is no 
consensus on how to express correctly the macroscopic 
constitutive relations solely out of microscopic consider- 
ations and under various boundary conditions or loading 
histories. This very basic interrogation was illustrated in 
a recent debate on how to understand the stress distri- 
bution below a sand pile and especially how to account 
for the dependence on preparation protocols [7]. A new 
mechanical approach was proposed based on the concept 
of "fragile matter" [5] and force chains propagation mod- 
elling [8]. But recent experiments have dismissed this 
approach and evidenced results more consistent with the 
traditional framework of general elasticity [9]. 

In this letter we obtain new and non-int\iitive results 
on stresses measurements at the bottom of a granular 
assembly confined in a rigid cylinder. Two kinds of ex- 
periments were performed. First, the mass at the bottom 
of the column is measured as a function of the granular 
material filling mass. Second, similar measurements are 
produced with an overweight on the top of the granular 
material. Similar experiments were previously made by 
Vanel et al. [10] who showed that contrarily to previous 
experimental reports, one could obtain quite reproducible 
data provided a good control of the packing fraction ho- 
mogeneity and a mobilization of all the friction forces at 
the walls in the upwards direction. The existence of an 
overshoot effect was clearly demonstrated then, but no 



interpretation nor systematic measurement were under- 
taken. Note that this overshoot effect is also an outcome 
of the "OSL theory" [8] based on a vault propagation 
modelling but an important feature of the theory was 
not evidenced experimentally [10] : the existence of stress 
oscillations. In the present report, we use the same exper- 
imental set-up as in [10] but with a noticeable change in 
the measurement procedure which we claim is a decisive 
improvement for the data interpretation. 

The grains are dry, noncohesive, and slightly polydis- 
perse (10%) glass beads of diameter d = 1.5 mm piled 
into a vertical steel cylinder. A brass overweight can 
be added on top of the granular material. Three differ- 
ent column radii are used: i? = 19 mm, 28 mm and 40 
mm. The static coefficient of friction fjig between the glass 
beads and the steel walls is measured using the sliding 
angle of a three-bead tripod. We found: jig = 0.25±0.02. 
Two other steel cylinders of radius i? = 19 mm are used: 
a rough one (/i^ = 0.285 ± 0.02), and a smoother one 
(/Ltg = 0.22 ± 0.02). The column is closed at the bottom 
by a movable piston (the diameter mismatch is 0.5 mm). 
A force probe of stiffness fc = 40000 N.m"'^ is located 
under the piston and can be moved at constant driving 
velocity V via a stepping motor (see Fig. 1 inset). 

The experiments were performed on three different pil- 
ings: a dense one (64.5 ±0.5%) obtained by rain-filling, a 
loose one (59.0 ±0.5%) obtained by using an inner cylin- 
der which is removed slowly after filling, and a medium 
one (62.0 ± 0.5%) obtained by vibrating a loose packing. 

For exxch data point, the same procedure is repeated: 
first the column is filled with a mass of beads M jm and 
prepared according to the desired average packing frac- 
tion u. Then, the vertical translational stage is moved 
down slowly by a step motor (resolution of 100 nm/stcp) 
and the mass on the gauge is monitored. From the mass 
versus vertical descent data, the apparent mass Ma at the 
bottom of the column is extracted. We use a very slow 
descent of the force probe at a velocity Vq = 1-5 A*m.s~^, 
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more than ten times smaller than the one used by Vanel 
et al. [10]. Then, a quasi-static situation is attained 
were all dynamical oscillations are suppressed. In the 
first stage of the descent, the granular material does not 
move (see Fig. 1 inset) . The stress probe yields a value : 
m{t) = mo — {k/g)Vot, where mo is the apparent mass at 
time t = just after the end of the filling procedure (note 
that the value of mo is highly irreproducible from one fill- 
ing to the other). When the friction at the walls reaches 
the Coulomb threshold, the material starts to slip and 
then gently stabilizes at the value V{t) ~ Vq. Then the 
apparent mass may decrease for a piston descent of an- 
other 50 /im or more (the precise value depends on the 
column size) before it reaches a stationary state. This 
later evolution could be due to mechanical restructuring 
of the piling caused by the slippage at the boundaries. 
This evolution to a plateau limit should not be mixed 
with a dilatancy effect, which would occur much later, 
for displacements of about a grain size and also would 
lead to an increase of the apparent mass [10]. We choose 
to make the apparent mass Ma measurements, just at 
the end of the linear relaxation regime (see arrow on Fig. 
1 inset). Therefore, this new procedure should ensure 
(i) that the static piling is at the threshold of slippage 
everywhere at the walls and (ii) that we measure the me- 
chanical properties of the initial piling in the absence of 
mechanical restructuring. These are major differences 
with the original procedure of Vanel et al. [10] who ex- 
tracted the apparent mass at rest, after relaxation from 
the plateau value. In cylindrical coordinates with ori- 
gin at the top surface and the cylinder axis being the 
z axis, the relation, at the slipping onset, between the 
shear stress arz and the horizontal stress <Trr at the walls 
is (7rzir=R, z) = iisarr{i'=R, z) [jig is the Coulomb static 
friction coefficient between the grains and the wall). 

The typical results obtained whcin the filling mass Mfm 
is varied are shown on Fig. 1. The apparent mass Ma 
saturates exponentially with Mjui. When an overweight 
equal to the saturation mass Msat is added on top of 
the granular material. Ala increases with Mfm, up to 
a maximum M^ax, then decreases slowly towards the 
saturation mass Msat- We also verified that if we add 
and remove slowly the overweight, we can recover the 
apparent mass corresponding to the original Janssen's 
saturation curve within the experimental uncertainties. 

The simple model which captures the physics of this 
saturation phenomenon was provided in 1895 by H.A. 
Janssen [11]. This model is based on the equilibrium of a 
granular slice taken at the onset of sliding and predicts a 
relation between the filling mass Mfm and the apparent 
mass at the bottom Ma of the form: 

Ma = M,,t(l - exp(-^)) with M,at = ^ (1) 




M, 



sat 



where p is the mass density of the gramilar material, and 
K is the Janssen parameter rendering the average hori- 
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FIG. 1: Apparent mass Ma vs. filling mass Mfm for a loose 
packing in the medium-rough 38 mm diameter column with- 
out an overweight (squares), and with a 80.5 g overweight 
(triangles). Inset: measured mass vs. translation stage dis- 
placement Azstage, for a loosc packing in the smooth 38 mm 
diameter column; the arrow indicates the measure point; the 
dotted line corresponds to the linear spring relaxation (slope 
-k/g). Sketch of the experimental display. 



zontal redirection of vertical stresses. From a mechanical 
point of view, a major simplification of this model comes 
from the assumption that the redirection parameter K 
would stay constant along the vertical direction. But on 
the other hand, it provides a clear and simple physical ex- 
planation for the existence of an effective screening length 
A = R/2Kiis above which the mass weighted at the bot- 
tom saturates. Many saturation profiles were measured 
for various packing fractions, columns sizes and friction 
coefficients. Janssen's model may look very elementary 
but when the apparent mass rescaled by the saturation 
mass is plotted clS cL function of the filling mass also 
rescaled by the saturation mass, we obtain a universal 
resettling of all data on a curve which is precisely the one 
predicted by Janssen : Ma/M^at = fiMjui/Mgai), with 
f{x) = 1— exp(— x) (see Fig. 2). We recall that each data 
point corresponds to an independent experiment. The 
results is a priori so surprising that we will try next to 
understand it in reference to a simple elastic modelling of 
the granular column. The differences in Msat measured 
for different densities means that the effective Janssen's 
coefficient K increases with the average packing fraction 
(Fig. 2 inset). In other words the efficiency of vertical 
stresses redirection is better for higher packing fractions 
and it is a very sensitive result {AK/K ~ bAv/v). We 
also obtain a Janssen's coefficient K independent of the 
radius R and of the friction coefficient jis (Fig. 2 inset). 
These results are to the best of our knowledge the first 
experimental validation of this 1895 model. 

For experiments with an overweight on the top, we also 
rescale the data for different column sizes (Fig. 3a) and 
also find a universal curve g{x) such that Ma/Mgai. = 
g{M fill/ Msat) which means that the data scale with 
as Msat does. Note that this overshoot goes beyond the 
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FIG. 2: Apparent mass Ma vs. filling mass Mfm, rcscalod by 
the saturation mass Msat, for a loose packing in the smooth 
38 mm diameter column (triangles), for a dense packing in 
the rough (open circles) and smooth (squares) 38 mm diame- 
ter columns, and for a dense packing in the 80 mm diameter 
column (open down triangles). The results are compared to 
Janssen's model prediction (line). Inset: Janssen's param- 
eter K extracted from experiments made in medium-rough 
columns of 3 diameters: 38 mm (squares), 56 mm (circles), 
80 mm (triangles), and in a smooth (open down triangles) 
and a rough (open diamonds) 38 mm diameter column, vs. 
packing fraction. 



possibilities of Janssen's model which would predict a flat 
curve: Ma/ Msat = 1- But, as wc cliangc density or fric- 
tion, we find no rescaling of the data. On Fig. 3b inset 
we display the maximum overshoot value Mmax rescaled 
by the saturation mass Mgat and show that its value in- 
creases with packing fraction and friction at the walls. 
Thus we can reformulate this result in a very surpris- 
ing and non-intuitive way: the more the weight of the 
grains is screened by the walls, the less the weight of the 
overload is screened\ Note that a E? scaling for the ap- 
parent masses (valid with and without overweight), rules 
out theoretical predictions [12] based on the q-model [13]. 
Furthermore, no stress oscillation is present contrary to 
the OSL theoretical predictions [10]. 

We now study in detail the predictions of isotropic ho- 
mogeneous elasticity. We consider an elastic medium, of 
Young modulus E and Poisson coefficient z/p, in a col- 
umn of radius R. Stresses and displacements can be cal- 
culated in the limit of high depths z under the assump- 
tion that they then should be independent of z. The 
boundary conditions we impose are the Coulomb condi- 
tion everywhere at the walls {(Jrz{'''=R) = iJ,s<7rr{r=R)), 
and infinitely rigid walls i.e. Ur{r = R) = 0. Using the 
stress-strain relation Ecij = (1 -|- VpjcJij — '-'pSijcrkk and 
internal equilibrium relation diaij = —pgj we find cr^r = 
(Tee = - i^p))cTzz with alf{r,z) = -f^. The 

asymptotic displacements arc ti^ (r, z) = — pgr"^ — 

^ 2^^iy PflRz + ^0 and z) = ug{r, z) = 0. Thus, we 
obtain a Janssen's like redirection phenomenon due to 



a Poisson's ratio effect with a local Janssen's parameter 
Kei = Urrir, z)/azz{r, z), being for large depths: 

Kei = //p/(l - lyp) (2) 

Thus, in the context of isotropic elasticity, the depen- 
dence of K on packing fraction can be interpreted as an 
increase of the Poisson's ratio Vp with packing fraction. 
Note a close derivation of stress distribution in an elastic 
column for R (this condition is not realized in our 
experiment) obtained by Evesque and de Gcnnes [14]. 

In order to get the whole stress saturation curve, finite 
element numerical simulations [15] were performed. The 
column is modelled as an isotropic elastic medium. We 
vary the friction at the walls pa, the Young modulus E 
and the Poisson coefficient Up. We imposed a rigid, ei- 
ther perfectly stick or perfectly slip bottom. We found no 
appreciable difference between these two previous cases. 
The condition arz = Ps<Trr is imposed everywhere at the 
walls. The cylinder is modelled as a steel elastic medium. 
As long as the Young modulus E of the elastic medium is 
less than 500 MPa, which is usually the case for granular 
media, we found no dependence of the results on E. We 
verified that in all the simulations performed, there is no 
traction in the elastic medium, so that this can actually 
be a model for a granular material. We observed indeed 
that for low friction coefficients i.e. ps < 0.5, the stress 
curve at the bottom of the column could not be sepa- 
rated from Janssen's universal curve. Note that the pres- 
ence of a bottom imply that the effective Janssen's pa- 
rameter K^i^ extracted from Janssen's scaling is higher 
than the Kej. Wc have in fact : K^f'^ = aK^i + b with 
a = 1.029 ± 0.002 and b = -0.008 ± 0.001. If we use the 
correspondence with the eff'ective Janssen's constants ob- 
tained experimentally this would yield an empirical rela- 
tion between Poisson ratio Vp and packing fraction v as: 
Up ~ 2.3(z/ — 0.41) with a precision of 5%. Note that the 
largest packing fraction u = 0.645 ± 0.005 would give a 
Poisson ratio Up = 0.54 ± 0.03 marginally larger than the 
limit value of 1/2. Note that i) nothing really insures 
that our piling preparation is strictly isotropic and ii) in 
spite of the careful procedure, we are never absolutely 
sure that all the friction forces at the wall are actually 
fully mobilized upwards. We also performed simulations 
with an overweight (which is taken perfectly stick or slip 
with no appreciable difference), using brass mechanical 
parameters. We found that the resulting curve is very 
similar to the experimental one, but the experimentally 
observed maxima are consistently 30 to 40 times larger 
(see Fig. 3b). This is what we call the giant overshoot 
effecf'l We find from numerical computations the same 
qualitative phenomenology as in the experiment i.e. the 
computed values of the overshoot Mmax rescaled by the 
saturation mass Msat-, increase both with the friction at 
the walls and with the Poisson coefficient (i.e. with the 
effective Janssen's parameter). Basically, wc recover in 
the elastic case, the same paradox as the experimental 
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FIG. 3: a: Apparent mass Ma vs. filling mass Mfm, rescaled 
by the saturation mass Msat, for loose packings in medium- 
rough columns of 3 diameters (38 mm (squares), 56 mm (cir- 
cles), 80 mm (triangles)) with an overweight equal to Msat', 
the dotted line is the hydrostatic curve, b: Simulation of the 
experiment of Fig. 3a, for an elastic medium characterized 
by the same saturation mass (Poisson coefficient Vp = 0.46) 
and the same friction at the walls (jis — 0.25). Inset: maxi- 
mum mass Mmax rescaled by saturation mass A4sat vs. static 
coefficient of friction at the walls fis , in experiments made on 
loose (squares) and dense (circles) packings in 38 mm diam- 
eter columns, and in simulations for elastic media of Poisson 
coefficients i/p = 0.45 (open squares) and Vp = 0.49 (open cir- 
cles) ; the left vertical scale is used for the experimental data, 
the right vertical scale is used for the simulation data. 

situation. We now understand it as a consequence of 
the boundary condition imposed experimentally by the 
overweight i.e. an almost constant displacement on the 
surface : Uz{r) = uo- When Vp or iig is decreased to 
zero (i.e when Msat increases), the parabolic part of the 
asymptotic vertical displacement becomes negligible and 
therefore, the displacement on the surface is close to the 
asymptotic value. This leads to a flat pressure profile. 

Understanding the amplitude of this overshoot effect 
is still an open question but we propose here a qualita- 
tive vision based on stress induced anisotropy. Small me- 
chanical restructuring due to displacements close to the 
overweight at the top could be capable to produce lo- 
cal anisotropy in a way that the vertical direction would 
becomes stiflFer due to the preferential increase in the 
number of granular contacts. A consequence of this 
anisotropy is that the overshoot amplitude can be in- 
creased. This is what we verified numerically using a 
model of orthotropic elasticity [15] with the stiff axis 



taken in the vertical direction. But importantly, we need 
also to assume that this structuring effect has a finite 
spatial extend since for large column heights the original 
saturation limit is recovered. 

In summary, we performed series of precise exper- 
iments on a granular material confined in cylindrical 
columns. We obtained for the first time the universal 
scaling predicted by Janssen in 1895 for the stress satu- 
ration curve when the filling mass is increased. We show 
that the screening parameter (the effective Janssen's con- 
stant) is an increasing function of the average packing 
fraction. By comparison with isotropic elastic theory we 
provide a one to one correspondence between the Pois- 
son's ratio and granular packing fraction. We study sys- 
tematically the overshoot effect, i.e. when a mass equal 
to the saturation mass is added to the top of the packing 
we have a non-monotonic variation of the apparent mass 
as a function of the filling mass. This effect is reproduced 
qualitatively by isotropic elasticity but of a scale 30 to 
40 times smaller. We suggest a qualitative interpretation 
of this effect based on stress induced anisotropy. There- 
fore along these lines, this giant overshoot effect can be 
seen as a strong test for any theoretical model describing 
the statics of granular assemblies under various loading 
histories and boundary conditions. 

We acknowledge discussions with P. Claudin and J. 
Socholar, and thank J. Lanuza for technical assistance. 
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